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Computationally  Efficient  Substrate  Noise  Coupling  Estimation  in  Lightly  Doped 

Silicon  Substrates 


1.  INTRODUCTION 


The  digital  revolution  has  driven  the  need  for  highly  integrated  analog,  RF  and 
digital  systems  to  co-exist  on  a  common  silicon  substrate.  While  this  integration  enhances 
speed  and  reduces  power  and  cost,  it  brings  with  it  the  issue  of  substrate  noise.  A 
simplified  representation  of  this  issue  is  shown  in  Figure  1.1.  The  digital  switching  noise 
couples  through  the  supply  lines  and  the  common  substrate  to  the  analog  and  RF  blocks. 
This  results  in  a  performance  degradation  of  the  analog  and  RF  blocks.  The  noise 
coupling  increases  with  an  increase  in  the  frequency  of  operation  and  also  with  a  decrease 
in  the  device  dimensions,  thus  making  substrate  noise  characterization  crucial  in  the 
design  cycle. 


Digital 


Analog 


Figure  1.1.  Substrate  noise  coupling  from  a  digital  circuit  to  an  analog  circuit. 
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The  noise  propagation  properties  of  the  substrate  are  strongly  dependent  on  its 
doping  levels.  Based  on  the  doping  concentration,  substrates  can  be  classified  into  two 
primary  types:  lightly  doped  (high  resistivity)  substrates  and  heavily  doped  (low 
resistivity)  substrates.  Latchup  issues  in  lightly  doped  substrates  resulted  in  an  increased 
popularity  of  epitaxial  or  heavily  doped  substrates.  With  technology  scaling  in  recent 
years,  the  trend  has  shifted  towards  lower  voltage  operations,  and  hence  latchup  has 
become  less  of  an  issue  [1].  In  addition,  an  epitaxial  substrate  proves  less  effective  in 
isolating  the  noisy  digital  blocks  from  the  sensitive  analog  cells  when  the  circuit  is 
densely  packed.  In  order  to  achieve  better  signal  integrity,  the  lightly  doped  substrate  is  a 
preferred  choice  for  today’s  ICs.  This  is  true  especially  for  applications  where  analog, 
RF  and  digital  circuitry  are  integrated  on  the  same  die  [2], 

Since  substrate  noise  is  capable  of  degrading  circuit  performance,  designers  are 
looking  for  better  ways  to  know  a-priori  the  noise  isolation  that  they  can  achieve  and 
techniques  to  improve  it.  Pre-layout  noise  analysis  at  the  schematic  or  cell  level  or  during 
floor  planning  can  be  more  beneficial  than  the  traditional  post-layout  analysis.  At  the  pre¬ 
layout  level,  a  quick,  yet  reliable  approximation  of  noise  coupling  between  the  noisy 
node  and  the  sensitive  node  is  required.  Numerical  methods  to  accurately  characterize  the 
substrate  are  readily  available.  Unfortunately,  numerical  methods  are  computationally 
expensive  and  thus  time  consuming.  This  makes  them  unsuitable  for  use  during  pre¬ 
layout.  Accuracy  can  be  traded-off  for  efficiency  during  the  pre-layout  noise  estimation 
stage  and  hence  there  is  a  thrust  to  find  quicker  but  less  accurate  solutions. 
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1.1.  Summary  of  Previous  Work  and  Contributions 

As  noted  earlier,  several  accurate  numerical  methods  characterize  the  lightly 
doped  substrate.  The  methods  can  be  classified  into  two  broad  categories:  finite 
difference  methods  [3],  [4]  and  boundary  element  methods  [5],  [6]. 

In  finite  difference  methods  (FDM)  the  entire  volume  of  the  substrate  and  the 
contacts  placed  on  it  are  divided  into  a  mesh  and  Poisson’s  equation  is  used  to  estimate 
the  potential  variation.  The  rigorous  3-D  meshing  makes  it  a  time  consuming  method.  In 
the  boundary  element  method,  the  surfaces  of  the  substrate  contacts  are  divided  into 
panels  to  solve  equations  for  the  potential  variation.  The  method  is  faster  than  FDM,  but 
is  still  not  a  good  choice  for  noise  estimation  at  the  pre-layout  stage. 

Empirical  models  are  alternatives  to  the  numerical  methods.  They  are  obtained  as 
a  result  of  a  curve  fitting  process  applied  to  simulation  or  measured  data  [7],  [8],  [9]. 
Such,  models  often  express  noise  coupling  in  terms  of  the  geometry  and  the  separation  of 
the  contacts.  Unfortunately,  empirical  methods  do  not  produce  generalized  models  and, 
hence,  cannot  be  applied  to  arbitrary  substrates.  In  addition,  several  assumptions  made 
while  developing  such  models  lead  to  accuracy  issues  when  applied  to  situations  where 
these  assumptions  are  violated.  However,  the  empirical  models  provide  a  quick 
estimation  of  noise  coupling  during  the  preliminary  stages  of  circuit  design. 

This  work  focuses  on  developing  a  simple  and  reliable  resistive  macromodel  that 
characterizes  the  lightly  doped  substrate  up  to  frequencies  of  about  a  few  gigahertz.  A 
brief  comparison  with  existing  approaches  and  a  summary  of  the  contributions  is: 
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•  In  [8],  the  effective  coupling  resistances  are  modeled  directly.  This  can 
lead  to  issues  when  extending  the  two-port  model  to  practical  multi¬ 
contact  situations.  This  work  adopts  a  Z-parameter  based  approach  to 
develop  a  general  model. 

•  The  model  described  in  [9]  has  unacceptable  errors  when  applied  to  very 
closely  spaced  contacts,  i.e.,  of  the  order  of  l~10pm.  The  model 
developed  in  this  thesis  provides  better  accuracy  for  these  situations. 

•  In  [7],  a  Z-parameter  based  model  has  been  developed.  The  model 
requires  several  process  constants  and  was  validated  for  a  few  contact 
geometries  and  separations.  The  following  work  attempts  to  develop  a 
compact  model  that  works  for  a  large  set  of  contact  geometries  and 
separations. 

1.2.  Thesis  Outline 

The  thesis  is  organized  as  follows:  Chapter  2  discusses  the  properties  of  lightly 
doped  substrates  and  the  reasons  for  choosing  open  circuit  parameters  to  model  the 
resistive  network.  In  Chapter  3,  the  model  development  for  the  cross-coupling  impedance 
and  the  self-impedance  is  described.  Chapter  4  explains  a  paneling  approach  and  the 
application  of  the  macromodel  in  conjunction  with  the  paneling  approach  to  account  for 
the  proximity  effects  of  neighboring  contacts.  Chapter  5  provides  comprehensive 
validation  of  the  macromodel  using  measurements  and  simulations.  The  macromodel 
with  paneling  is  applied  to  a  multi-contact  example  to  prove  the  usefulness  of  this  work. 
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Finally  Chapter  6  concludes  the  current  work,  and  discusses  the  scope  for  future 
investigations. 
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2.  PROPERTIES  OF  THE  LIGHTLY  DOPED  SUBSTRATE  AND  THE 
SUBSTRATE  COUPLING  NETWORK 

In  this  chapter  a  discussion  of  the  conductive  properties  of  a  non-uniform  lightly 
doped  substrate  is  provided.  This  is  followed  by  a  brief  explanation  for  representing  the 
lightly  doped  substrate  with  an  equivalent  network.  Open-circuit  parameters  are  used  to 
characterize  the  substrate  and  hence  the  noise  coupling.  Measuring  open-circuit 
parameters  requires  a  reference.  It  is  assumed  that  the  backplane  of  the  substrate  can  be 
chosen  as  the  reference.  The  conditions  under  which  this  assumption  is  valid  are 
discussed. 

2.1.  Conductive  properties  of  lightly  doped  substrates 

The  lightly  doped  substrate  is  a  high  resistive  substrate  with  bulk  resistivities  of 
the  order  of  10~50Q-cm.  Noise  isolation  can  be  improved  by  simply  increasing  the 
distance  between  the  noisy  and  sensitive  blocks  [6],  However  the  high  resistivity  of  such 
substrates  is  often  decreased  because  of  p+  surface  implants  (for  p-type  substrates)  such 
as  channel-stop  implants.  Channel-stop  implants  are  usually  added  to  prevent  the 
formation  of  parasitic  paths  [10],  A  parasitic  conductive  layer  may  be  formed  when  a 
high  voltage  is  applied  on  a  uniformly  doped  silicon  surface  (of  p-type)  without  such 
implants.  The  conductivity  of  these  implants  is  about  two  orders  of  magnitude  more  than 
the  bulk.  The  implant  depth  is  typically  1pm.  The  higher  conductivity  aids  in  increasing 
the  threshold  voltage  of  the  surface  and  prevents  a  parasitic  conduction  path.  However, 
the  channel-stop  implant  makes  the  substrate  more  conductive  at  the  surface  and, 
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therefore,  analog  or  RF  circuits  are  not  isolated  to  the  fullest  extent  from  the  digital 
blocks. 

To  illustrate  the  above  effect,  the  equivalent  resistances  for  two  contacts  in  a 
uniformly  doped  substrate  and  a  substrate  with  a  channel-stop  implant  are  compared.  The 
equivalent  resistance  is  measured  by  injecting  a  unit  current  through  one  contact  and 
measuring  the  potential  on  the  other.  The  cross  sections  of  the  substrate  profiles  are 
shown  in  Figure  2.1.  The  resistances  were  extracted  from  a  Green’s  function-based 
solver,  EPIC  [5].  An  overview  of  the  tool  capabilities  and  assumptions  are  given  in 
APPENDIX  A.  The  plot  of  the  equivalent  coupling  resistance  between  the  two  contacts 
as  a  function  of  separation  is  presented  in  Figure  2.2.  The  coupling  resistance  is  an  order 
of  magnitude  lower  in  the  substrate  with  channel-stop  implants.  All  the  substrates 
reported  in  this  work  are  assumed  to  have  a  channel-stop  implant. 


p-type 

substrate 


lpm 


250pm 


(a) 


p+  channel-stop 
implant 


p-type  substrate 


Figure  2.1.  Cross  sections  of  lightly  doped  substrates,  (a)  Uniform  lightly  doped 
substrate,  (b)  Lightly  doped  substrate  with  channel-stop  implant. 
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Figure  2.2.  Comparison  of  the  equivalent  coupling  resistances  for  a  lightly  doped 
substrate  with  and  without  a  channel-stop  implant. 


2.2.  Equivalent  substrate  network  for  lightly  doped  substrates 

For  frequencies  less  than  a  few  gigahertz,  the  coupling  between  two  contacts  in  a 
substrate  can  be  represented  equivalently  by  a  resistive  7t-model  as  shown  in  Figure  2.3 
[7],  [8],  [9].  Rn  represents  the  cross-coupling  resistance  and  Rn  or  R22  characterize  the 
mutual  resistance  between  the  contact  and  the  backplane. 


Cl  C2 


Figure  2.3.  A  71-resistive  model. 
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The  open  circuit  or  Z-parameters  are  used  to  characterize  the  substrate  network. 
To  illustrate  the  open  circuit  parameters,  consider  in  general  a  set  of  contacts  on  the 
substrate.  If  V  is  the  vector  of  voltages  applied  on  contacts,  I  is  the  vector  of  currents  that 
is  conducted  by  the  substrate  and  Z  is  the  impedance  matrix,  then  the  current-voltage 
relationship  is  given  by 

[V]  =  [Z]*[I]  (2.1) 

For  a  simple  case  of  two  contacts,  the  above  relation  can  be  written  as: 


where  Z„  is  the  self-impedance,  or  the  potential  at  contact  Ci  when  a  unit  current  is 
injected  into  it  while  the  other  contact  Cj  (j^i)  is  kept  floating.  Zy  (f£j)  is  the  cross¬ 
coupling  impedance,  or  the  potential  on  a  floating  contact  Ci,  when  a  unit  current  is 
injected  into  the  other  contact  Cj.  The  impedances  are  given  by  the  following 


expressions: 


(2.3) 


By  the  principle  of  reciprocity,  Zy=Zji.  The  Z-parameters  are  a  function  of  the  sizes  and 
separation  of  the  two  contacts  for  a  given  substrate.  This  fact  is  used  to  model  these 
parameters  in  this  work.  The  Z-parameters  can  be  related  to  the  resistive  network  shown 
in  Figure  2.3  by  the  following  relations: 
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Open  circuit  parameters  are  chosen  to  model  the  substrate  because  they  are 
largely  independent  of  the  surrounding  contacts.  This  follows  from  the  fact  that  the 
surrounding  contacts  are  open  circuited  or  left  floating  while  measuring  the  self¬ 
impedance  for  a  contact  (C 1  or  C2)  and  the  cross-coupling  impedance  between  the  two 
contacts  (C 1  and  C2).  The  two-contact  problem  can  be  easily  extended  to  multi-contacts 
by  individually  computing  the  Z  matrix  for  every  pair  of  contacts.  This  property  is  not 
exhibited  by  the  resistance  matrix  or  the  admittance  matrix,  which  makes  them  unsuitable 
to  characterize  the  substrate  network  [11].  Finally,  to  extract  the  resistance  network  for  a 
multi-contact  case,  the  Z-matrix  can  be  inverted  to  obtain  the  admittance  matrix  and  the 
resistances  can  be  calculated  from  this  admittance  matrix  (Section  4.2). 

Characterization  of  a  lightly  doped  substrate  is  a  more  involved  problem 
compared  to  characterizing  a  heavily  doped  substrate.  After  a  certain  threshold 
separation,  a  heavily  doped  substrate  conducts  mainly  through  its  bulk  and  not  through 
the  surface.  Beyond  this  threshold  separation,  the  cross-coupling  resistance  between  any 
two  contacts  becomes  too  large  and  the  cross-coupling  impedance  becomes  negligible 
(from  Equation  (2.4)).  Hence,  there  is  a  specific  range  of  separations  for  which 
characterization  is  required.  Due  to  the  presence  of  a  channel-stop  implant,  conduction 
takes  place  through  the  surface  in  a  lightly  doped  substrate  even  for  large  separations. 
This  indicates  that  the  coupling  resistance  is  not  large  enough  at  far  separations  and, 
hence,  the  cross-coupling  impedance  cannot  be  neglected  [6]. 
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2.3.  Backplane  properties  of  a  lightly  doped  substrate 

Open  circuit  parameters  require  a  reference  for  their  measurements.  For  a 
substrate  network,  the  backplane  is  often  chosen  to  be  the  reference.  Consequently  the 
backplane  would  have  to  be  assumed  equipotential.  The  backplane  is  usually  not 
considered  equipotential  for  a  lightly  doped  substrate  [12].  However,  the  backplane  can 
be  approximated  as  an  equipotential  surface  for  certain  conditions  as  explained  in  this 
section. 

If  the  backplane  of  a  substrate  is  floating  the  equivalent  network  for  a  two- 
contact  case  consists  simply  of  an  effective  cross-coupling  resistance  (Ri2f).  This  cross¬ 
coupling  resistance  (Ri2f)  can  be  extracted  from  the  grounded  backplane  network 
provided  the  backplane  is  an  equipotential  surface  (Equation  (2.5))  as  shown  by  Figure 


2.4. 


RUf  -  Rn  II  (Rn  +R22) 
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Figure  2.4.  Floating  backplane  substrate  network  derived  from  a  71-resistive  network. 
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If  the  relation  in  Equation  (2.5)  for  two  contacts  is  valid  for  the  lightly  doped 
substrate  under  consideration  (Figure  2.2(b))  then  the  backplane  can  be  approximated  as 
a  single  node.  First  the  above  equation  is  validated  for  two-contact  cases.  Then  multi 
contact  examples  are  considered.  For  the  latter  case,  SPICE  is  used  to  obtain  the  effective 
floating  cross-coupling  resistances  from  the  grounded  backplane  network  and  compared 
with  the  actual  floating  backplane  substrate  resistances.  All  the  above  validations  were 
conducted  using  EPIC  simulations  using  an  option  for  grounding  or  floating  the  substrate 
backplane. 

The  results  obtained  for  one  of  the  two-contact  cases  is  shown  in  Figure  2.5. 


Figure  2.5.  Comparison  of  the  effective  cross-coupling  resistances  between  two  contacts 
from  a  floating  backplane  setup  and  those  derived  from  a  grounded  backplane. 

The  comparison  between  the  backplane  floating  substrate  resistances,  and  the 
resistances  derived  from  the  backplane  grounded  substrates  for  the  multi-contact  case  is 
made  in  Table  2.1.  The  contact  layout  for  this  example  is  given  in  Figure  2.6.  In  both  the 
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experiments  the  relative  errors  for  all  cases  are  less  than  1%  when  the  effective  coupling 
resistances  are  compared. 
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Figure  2.6.  Layout  of  5  non-identical  contacts  with  dimensions  (1)  5pm  x  5pm,  (2)  15 pm 
x  22.5pm,  (3)  5pm  x  15pm,  and (4)  2.5pm  x  12pm,  and  (5)  3pm  x  8pm. 


Table  2.1.  Comparison  of  the  effective  cross-coupling  resistances  from  a  floating 
backplane  setup  and  those  derived  from  a  grounded  backplane  for  a  multi-contact  case. 


Rij 

From  floating  backplane  (Q) 

Extracted  from  grounded  backplane  (Q) 

Rl2 

277.6 

277.56 

R23 

369.56 

369.42 

R34 

397.03 

396.96 

R45 

398.05 

397.99 

14 


Another  example  of  a  multi-contact  case  with  large  contacts  placed  at  various 
locations  of  a  1mm  x  lnnn  die  is  presented  in  Table  B-l  of  Appendix  B.  The  maximum 
relative  errors  are  about  16%.  This  serves  as  yet  another  proof  for  the  claim  that  the 
backplane  of  this  substrate  can  be  considered  equipotential  with  a  reasonable  accuracy.  In 
addition,  the  same  multi-contact  example  is  considered  to  validate  the  equipotential 
backplane  assumption  for  a  resistive  ground  connection  and  is  also  presented  in 
Appendix  B,  following  the  previous  experiment.  From  the  analyses  of  two-contact  and 
multi-contact  examples,  it  can  be  concluded  that  the  backplane  of  this  lightly  doped 
substrate  can  be  approximated  as  an  equipotential  surface  with  a  reasonable  accuracy. 

The  substrate  thickness  plays  an  important  role  in  the  equipotential 
approximation.  In  order  to  analyze  the  backplane  characteristics  for  lightly  doped 
substrates,  tests  were  conducted  with  varying  bulk  thicknesses  from  10pm  to  250pm 
using  the  MEDICI  device  simulator  [13].  A  die  width  of  240pm  was  assumed,  as  only  a 
limited  number  of  grid  points  could  be  used  in  MEDICI.  Three  contacts  were  considered 
on  the  substrate’s  surface.  A  potential  of  IV  was  applied  on  one  contact  placed  in  the 
middle  of  the  die  while  grounding  the  other  two  contacts  (Figure  2.7).  The  potential 
variation  from  one  edge  of  the  backplane  to  the  other  was  measured  for  all  substrates. 
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Figure  2.7.  Simulation  setup  for  measuring  the  backplane  potential  variation. 

Assume  that  for  any  applied  voltage  x,  Ax  is  the  potential  variation  across  the 
backplane.  It  was  observed  that  there  was  a  linear  relationship  between  the  voltage 
applied  and  the  voltage  at  the  backplane  (Ax=kx).  Hence  the  reference  for  measuring  the 
potential  variation  across  the  backplane  was  taken  to  be  the  applied  voltage.  The 
backplane  potential  variation  relative  to  the  applied  voltage  for  different  substrate 
thicknesses  is  summarized  in  Table  2.2.  The  potential  variation  across  the  backplane  plots 
is  shown  for  a  250pm  thick  substrate  in  Figure  2.8  and  for  a  10pm  thick  substrate  in 
Figure  2.9.  The  potential  variation  for  a  10pm  thick  substrate  is  around  90%  whereas  it  is 
approximately  10%  for  the  250pm  thick  substrate,  with  respect  to  the  applied  voltage. 
From  the  above  observations  it  is  clear  that  the  potential  across  the  backplane  varies  by 
significant  amounts  only  for  substrates  with  thicknesses  less  than  200pm. 
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Table  2.2.  Potential  variation  across  the  backplane  for  different  substrate 
thicknesses  relative  to  the  applied  voltage. 


Substrate  thickness  (pm) 

Relative  backplane  potential  variation 

10 

90% 

100 

50-60% 

200 

18-20% 

250 

1-10% 

The  reason  behind  such  a  variation  can  be  explained  as  follows.  For  thicker 
substrates  an  electric  charge  placed  on  the  substrate’s  surface  views  the  entire  backplane 
as  a  point  object  and  hence  the  variation  in  potential  across  the  whole  backplane  is 
minimized.  After  observing  the  potential  variation  for  the  various  substrate  thicknesses  it 
was  concluded  that  the  assumption  of  an  equipotential  backplane  is  valid  for  substrates 
with  thicknesses  greater  than  200pm.  To  further  demonstrate  that  the  backplane  of 
thinner  substrates  cannot  be  considered  equipotential,  a  100pm  thick  substrate  with  five 
large  contacts  placed  on  it  is  considered  in  Appendix  B.  The  effective  floating  cross¬ 
coupling  resistances  derived  from  a  grounded  backplane  and  the  actual  floating 
backplane  resistances  are  compared  in  Table  B-2.  EPIC  was  used  for  simulating  the 
networks.  The  maximum  relative  errors  are  of  the  order  of  40%  when  this  backplane  was 
considered  as  an  equipotential  surface. 
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Figure  2.8.  Potential  variation  across  the  backplane  for  a  250pm  thick  substrate. 


Figure  2.9.  Potential  variation  across  the  backplane  for  a  10pm  thick  substrate. 
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3.  MACROMODEL  FOR  LIGHTLY  DOPED  SILICON  SUBSTRATES 

In  Chapter  2,  a  brief  overview  of  the  conductive  properties  of  a  lightly  doped 
substrate  was  provided.  The  assumption  that  the  backplane  of  the  substrate  can  be 
considered  to  be  an  equipotential  surface  was  validated.  It  was  also  noted  that  the  Z- 
parameters  of  the  substrate  depend  on  the  geometry  and  separation  of  the  contacts.  Given 
this  background,  the  goal  of  this  work  is  to  develop  a  model  that  expresses  the  Z- 
parameters  as  a  function  of  geometry  and  separation  of  the  contacts. 

This  chapter  focuses  on  the  development  of  a  model  for  the  Z-parameters  of  a 
lightly  doped  substrate.  The  model  is  initially  developed  for  a  simple  case  of  two  contacts 
for  a  typical  0.25pm  process.  EPIC  is  used  to  simulate  the  Z-parameters  for  several 
contact  cases  and  provides  the  basis  for  model  development  and  validation. 

3.1.  Existing  models 

There  are  two  main  parametric  approaches  to  model  the  noise  coupling  for  lightly 
doped  substrates.  These  approaches  differ  in  the  manner  in  which  the  geometry  and 
separation  of  contacts  are  quantified. 

In  [7],  the  relative  position  of  contacts  is  described  separately  along  the  X  and  Y 
axes.  A  drawback  of  this  approach  is  that  several  process  constants  are  needed  to  model 
the  Z-parameters.  In  addition,  effects  such  as  proximity  of  neighboring  contacts  are  not 
accounted  when  applying  the  model  to  multiple  contacts.  As  a  result,  this  model  may  not 
be  adequate  for  large  examples. 
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The  model  developed  in  [9]  employs  the  geometric  mean  distance  (GMD)  to 
measure  the  contact  separation.  The  GMD  based  approach  eliminates  the  use  of  separate 
parameters  to  define  the  contact  positions  along  the  X  and  Y  axes.  The  model  is  very 
efficient,  especially  for  arbitrarily  placed  contacts.  A  single  parameter  (GMD)  describes 
the  relative  positions  of  the  contacts  and  as  a  result  the  number  of  process  parameters  is 
minimized.  However,  the  geometric  mean  distance  suffers  from  certain  limitations 
leading  to  a  significant  loss  of  accuracy  in  practical  situations,  as  will  be  discussed  in 
Section  3.3. 

3.2.  Assumptions  for  the  macromodel  development 

In  Chapter  2,  preliminary  simulations  showed  that  the  effective  circuit  parameters 
for  a  floating  backplane  substrate  can  be  derived  from  the  Z-parameters  of  a  grounded 
backplane  substrate.  It  was  also  validated  that  the  backplane  of  substrates  with 
thicknesses  >  200pm  can  be  considered  equipotential.  A  typical  high  resistive  substrate 
closely  depicting  a  0.25pm  process  is  shown  in  Figure  3.1.  Assuming  that  the  substrate  is 
grounded,  EPIC  was  used  to  extract  the  Z-parameters  for  a  large  set  of  two  contact 
geometries  (115  test  cases  listed  in  Table  C-l  of  Appendix  C).  The  contact  dimensions  in 
the  test  cases  range  from  lpm  to  200pm.  A  large  variation  in  the  data  set  was  chosen  to 
ensure  that  the  model  will  work  for  a  large  set  of  geometries.  The  die  size  is  assumed  to 
be  lmmxlmm.  The  test  cases  consider  different  contact  geometries  that  include  squares 
and  rectangles  placed  in  different  orientations.  The  types  of  contacts  and  orientations 
considered  are  described  in  Appendix  C. 


20 
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250pm 


p+  channel-stop  implant 


p-type  substrate 


Figure  3.1.  Cross  section  of  a  lightly  doped  substrate  for  a  typical  0.25pm  process. 


The  coupling  trends  between  two  contacts  in  the  substrate  are  analyzed  by 
systematically  placing  the  contacts  in  several  locations  on  the  die  by  either  increasing 
their  edge-to-edge  or  center-to-center  separations.  The  definitions  of  these  separations 
between  a  pair  of  contacts  are  shown  in  Figure  3.2. 


fcj 

r 

edge-to-edge  separation 


center-to-center  separation 


(a)  (b) 

Figure  3.2.  Convention  ofXee  and  Ycc  separation  between  contacts,  (a)  Edge-to-edge 
separation  or  Xee  separation,  (b)  Center-to-center  separation  or  Ycc  separation. 

The  distance  Xee  denotes  the  horizontal  separation  between  the  closest  edges  of 

the  contacts,  i.e.,  the  edge-to-edge  separation.  The  distance  Ycc  represents  the  vertical 

separation  between  the  centers  of  the  contacts,  i.e.,  center-to-center  separation. 
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3.3.  Geometric  Mean  Distance  (GMD) 

The  geometric  mean  distance  (GMD)  was  introduced  in  [8]  to  measure  the 
effective  separation  between  contacts.  In  [8],  the  GMD  was  used  to  model  the  cross¬ 
coupling  resistance  between  the  contacts.  A  similar  approach  to  model  the  cross-coupling 
impedance  or  Zn  is  followed  in  this  thesis.  One  significant  advantage  of  the  GMD  is  that 
it  simultaneously  incorporates  the  effects  of  separation  between  contacts,  their  sizes  and 
orientation.  The  GMD  eliminates  the  need  to  individually  consider  the  effect  of  these 
parameters. 

The  definition  of  geometric  mean  distance  between  two  contacts  is  given  by: 


In  d 


GMD 


In  rdw,dLdw2dL 

JL2  J W2  J Li  J W,  ii// 

WxLxW2L2 


(3.1) 


where  Jgmd  is  the  geometric  mean  distance  and  r  is  the  Euclidean  distance  between 
arbitrary  points  (wi,  f)  and  (w2,  I2)  on  the  contacts  Cl  and  C2,  respectively  as  shown  in 
Figure  3.3  [14].  Wi,  Li,  W2,  L2  are  the  widths  and  lengths  of  contacts  Cl  and  C2, 
respectively.  The  actual  expression  can  be  simplified  [8],  and  is  of  the  fonn: 


gmd  =  (dld2 . dXb )16  (3.2) 

where  d,  is  the  Euclidean  distance  between  the  midpoints  of  the  edges  of  the  contacts  (8 
midpoints  yielding  16  distances). 
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Figure  3.3.  Geometric  mean  distance  between  two  contacts  using  midpoints. 


3.3.1.  Limitations  of  the  GMD 

Experiments  similar  to  the  ones  performed  in  [15]  have  been  conducted  for  a 
given  pair  of  contacts  in  a  lightly  doped  substrate.  First  the  Zi2  dependence  on  GMD  is 
checked  to  see  if  it  is  a  unique  function  of  the  GMD.  Then  conditions  under  which  ZJ2  is 
a  unique  function  of  the  GMD  are  determined. 

In  order  to  determine  how  Z12  is  influenced  by  the  GMD,  the  first  set  of 
simulations  consists  of  two  experiments: 

•  Experiment  1:  Analyze  the  relationship  between  Z12  and  GMD  for  a  two 
contact  case  by  varying  Xee,  keeping  Ycc= 0.  In  other  words,  the  centers  of 
the  contacts  are  aligned  (Figure  3.2(a)). 

•  Experiment  2:  Analyze  the  relationship  between  Z12  and  GMD  for  a  two 
contact  case  by  varying  Ycc  between  the  contacts  and  keeping  Xee  fixed  at 
2pm  (Figure  3.2(b)). 

The  results  from  these  experiments  are  shown  in  Figure  3.4.  Assume  that  the 
cross-coupling  impedance  is  a  function  of  the  GMD  alone.  Then  similar  GMD  values 
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from  these  experiments  must  generate  identical  Zn  values.  According  to  Figure  3.4,  at 
large  values  of  GMD  Zn  values  obtained  from  both  the  experiments  coincide.  However, 
at  closer  separations  (indicated  by  small  values  of  GMD),  Zn  values  from  both 
experiments  do  not  match. 


Figure  3.4.  Plot  of  Zn  values  from  EPIC,  as  a  function  of  GMD  for  2  contacts  (1pm  x 

1 10pm  and  5pm  x  27pm). 


Figure  3.5  shows  one  such  instance  where  identical  GMD  values  do  not  result  in  similar 
Zn  values.  The  GMD  is  identical  in  both  cases,  but  the  Zn  is  lower  by  17%  in  the  first 
case  (Experiment  1)  compared  to  the  second  case  (Experiment  2).  A  possible  explanation 
for  this  discrepancy  is  as  follows:  in  Case  2,  the  contacts  tend  to  remain  in  proximity, 
even  though  Ycc  is  changed.  However,  this  is  not  the  situation  in  Case  1.  The  proximity 
effects  influence  Zn,  in  addition  to  the  GMD.  As  a  result,  Z n  tends  to  be  higher.  Hence, 
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for  contacts  placed  at  very  close  separations,  GMD  alone  cannot  model  the  behavior  of 
Zl2- 


lpm 


1 10pm 


A^=15pm 


27pm 


5  pm 

Zi2=195.6Q 

GMD=33.7pm 


lpm 


1 10pm 


5  pm 


2jm 


27pm 

. A 


ycc=26fim 


Z12=230.0in 

GMD=33.7pm 


Case  1  Case  2 

Figure  3.5.  Example  to  show  identical  GMD  representing  different  cross-coupling 

impedance  values. 


In  the  above  analysis,  it  is  seen  that  there  exist  certain  cases  when  Z12  is  not 
dependent  on  the  GMD  alone.  Therefore,  specific  regions  where  Z12  is  not  dependent 
solely  on  GMD  need  to  be  identified,  for  a  given  process. 

In  Experiment  2,  theXee  separation  is  varied  at  regular  intervals,  say,  2pm,  10pm, 
20pm,  50pm.  For  each  Xee  separation,  the  Ycc  distance  is  varied  and  the  cross-coupling 
impedances  are  extracted  from  EPIC.  A  plot  of  Z[2  as  a  function  of  GMD  is  obtained  for 
each  Xee  separation  and  is  shown  in  Figure  3.6.  The  Z12  vs.  GMD  obtained  from 
Experiment  1  is  plotted  on  the  same  figure.  If  the  GMD  uniquely  defines  Z12  at  all  Xee  and 
Ycc  separations,  then  all  the  plots  would  coincide.  However,  this  is  not  the  case  and  the 
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relative  errors  for  Zn  (compared  to  the  Zn  values  I'rom  Xee  separation  with  7CC=0)  reduce 
from  17%  to  close  to  0%  as  Xee  increases.  It  is  seen  that  Zn  is  a  unique  function  of  the 
GMD  for  separations  where  Xee>  20pm,  for  the  given  substrate  and  contacts. 


Figure  3.6.  Z n  as  a  function  of  GMD  for  different  Xee  separations  when  Ycc  is  varied. 


When  the  above  experiments  were  repeated  for  contacts  with  different 
geometries,  it  was  observed  that  the  errors  were  more  pronounced  as  the  length  of  the 
contacts  increased  for  relatively  small  Xee.  Clearly,  the  geometric  mean  distance  is  not 
adequate  for  correctly  modeling  the  cross-coupling  impedances  at  very  close  separations. 
The  GMD  should  be  used  with  a  correction  parameter  to  model  all  regions  accurately. 
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3.4.  Cross-coupling  impedance  (Zi2)  model 

With  the  understanding  gained  from  the  previous  section,  the  modeling  of  Z12  for 
a  simple  case  of  two  contacts  is  discussed.  The  model  development  is  decomposed  into 
three  steps.  They  are  as  follows: 

i.  Assume  Ycc  =  0,  vary  Xee,  and  model  Z\ 2  as  a  function  of  GMD.  In  other  words, 
model  the  impedances  when  contact  centers  are  aligned. 

ii.  Next,  assume  Ycc  ±  0  and  use  the  above  GMD  model  to  characterize  the  cross¬ 
coupling  impedance.  In  other  words,  model  the  Z12  when  contact  centers  are  mis¬ 
aligned 

iii.  Use  a  correction  parameter  to  model  the  regions  where  GMD  does  not 
characterize  the  cross-coupling  impedance  between  contacts  correctly. 


3.4.1.  Z12  model  for  contacts  with  centers  aligned 

The  Z12  is  extracted  by  placing  the  contacts  in  the  center  of  the  die  and  Xee  is 
varied  from  1pm  to  600pm  keeping  Ycc= 0  (Figure  3.2(a)).  The  range  1pm  to  600pm  is 
chosen  since  the  largest  contacts  in  the  set  are  squares  with  a  dimension  of  200pm.  The 
squares  can  be  separated  by  a  maximum  of  600pm  in  a  lnnnx  1mm  die. 

Using  the  extracted  EPIC  results,  an  exponentially  decreasing  function  of  the 
geometric  mean  distance  is  used  to  model  the  mutual  coupling.  An  expression  for 
Z12  (centers  aligned)  that  follows  the  EPIC  trend  is  as  follows: 


^i2(c-c)  a\e 


£3 gmd 


(3.3) 
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where  ai  is  a  function  of  the  geometry  of  the  two  contacts  and  ks=0.  /9(uin)"((l  51  is  a 
process  constant.  The  curve  fitting  procedure  is  outlined  in  APPENDIX  D.  The  relative 
errors  between  the  EPIC  simulations  and  model  results  were  found  to  be  within  ±10%  for 
a  large  set  of  geometries  and  separation.  The  comparison  between  EPIC  simulations  and 
the  GMD  model  equation  for  four  different  cases  is  shown  in  Figure  3.7. 


(c) 


(d) 


Figure  3.7.  Comparison  of  the  GMD  model  for  Zn  and  EPIC  simulation  results.  The 
contact  dimensions  are  (a)  3pm  x  2.4pm  and  2pm  x  1pm,  (b)  5pm  x  40pm  and  5pm  x 
40pm,  (c)  100pm  x  100pm  and  100pm  x  100pm,  and  (d)  60pm  x  2pm  and  45pm  x 

lpm. 


Next,  the  parameter  ai  has  to  be  modeled  as  a  function  of  geometry.  The 
parameter  ai  cannot  be  modeled  as  the  merged  self-impedance  as  suggested  in  [9]  based 
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on  the  trends  observed  for  different  contact  geometries  (Figure  3.8(a)).  The  parameter  ai 
is  modeled  as  a  linear  function  of  the  sum  of  widths  of  the  two  contacts  as  given  below. 

a]  =  ( sum _of  _  widths)  +  k2  (3 .4) 

where  k}  and  k2  are  process  constants  and  sum_of_widths  is  the  sum  of  wj  and  w?  (Figure 
3.3).  This  equation  fits  the  aj  values  extracted  from  the  model  with  a  relative  error  of 
±6%  as  shown  in  Figure  3.8(b). 


Test  Cases 


(a) 


Figure  3.8.  a;  as  a  function  of  geometry,  (a)  Comparison  of  the  model  with  extracted 
values  from  the  GMD  model,  (b)  Relative  errors  for  the  aj  model  compared  to  the 
extracted  values  from  the  GMD  model. 
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Since  the  GMD  encapsulates  the  dependence  of  the  cross-coupling  impedance  on 
the  size  and  orientation  of  contacts,  a2  is  a  simpler  function.  Given  this  a2  model,  the 
cross-coupling  impedance  of  two  contacts  when  their  centers  are  aligned  ( Zn(c-c, ))  can  be 
represented  by  the  following  equation. 

Zi2(c-o  =  iki(wi  +  w2)  +  k2]e~k3^  (3.5) 

where  kj,  k2  and  k3  are  process  constants  and  w/  and  w2  are  the  widths  of  the  two  contacts 
(Figure  3.3).  To  summarize,  the  parameters  of  this  model  are  extracted  by  a  two-step 
curve  fitting  approach: 

1.  In  the  first  step,  Zl2  is  modeled  as  a  function  of  GMD.  Values  of  a2  are  obtained 
for  different  geometries  from  curve  fitting.  k2  is  extracted  as  a  process  constant. 

2.  In  the  second  step,  a2  is  modeled  as  a  function  of  the  contact  geometries.  ki  and  k2 
are  extracted  as  process  constants,  from  the  curve  fitting  process. 

Equation  (3.5)  can  be  directly  used  to  fit  the  Zn  simulation  data  instead  of  the  two- 
step  curve  fitting  and  reduce  the  errors.  From  such  a  direct  curve  fitting  approach  the 
values  of  the  process  constants  are: 

ki  =0.2322  D  (pm)'1  k2  = 638  Q  k3  =0.195  (pm)'05 

The  relative  errors  between  the  model  and  EPIC  simulations  are  within  ±10%  for  a  large 
range  of  geometries  and  separations.  The  relative  errors  increase  more  than  10%  for  very 
small  contacts  (which  are  below  a  dimension  of  5  pm)  at  separations  less  than  5  pm.  The 
model  is  also  limited  to  contact  geometries  with  length  ratios  (L1/L2  or  L2/L1)  less  than 
100.  The  comparison  of  the  model  equation  with  EPIC  values  for  another  four  cases  is 
shown  in  Figure  3.9,  where  the  entire  model  (Equation  (3.5))  has  been  employed. 
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Figure  3.9.  Comparison  of  the  GMD  model  (Equation  (3.5))  with  EPIC  simulations.  The 
contact  dimensions  are  (a)  1.2pmx  1.2pm  and  1.2pmx  1.2pm,  (b)  25pmx7pm  and 
16pmx20pm,  (c)  50pmx50pm  and  50pmx50pm,  and  (d)  125pmxlpm  and  125pmxlpm. 


A  compact  equation  with  just  3  process  constants  has  been  developed  to  model 
the  cross-coupling  impedance  between  contacts  that  have  aligned  centers.  The  next  step 
is  to  model  the  cross-coupling  impedance  for  contacts  where  the  centers  are  not  aligned. 


3.4.2.  Z12  model  for  contacts  with  centers  not  aligned 

In  practical  layouts,  the  contacts  or  ports  are  placed  arbitrarily,  often  with  their 
centers  not  aligned.  The  model  developed  in  the  previous  section  is  reasonably  accurate 
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when  the  centers  are  aligned.  The  subsequent  step  is  to  evaluate  the  performance  of  the 
model  equation  (Equation  (3.5))  by  applying  it  to  contact  cases  where  the  centers  are 
misaligned.  Impedances  were  extracted  from  EPIC  using  the  same  geometric  data  set,  but 
with  Ycc  varied  from  Opm  to  400pm  and  Xee  held  constant  (Figure  3.2(b)).  The  process  is 
repeated  for  different  Xee  (4pm,  6pm,  10pm  and  50pm)  to  cover  contact  cases  from  near 
field  to  far  field. 


(a)  (b) 


(c) 


(d) 


Figure  3.10.  Behavior  of  mutual  coupling  with  Ycc  separation.  The  contact  dimensions  are 
(a)  5pm  x  5pm  and  5pm  x  5pm,  (b)  20pm  x  1pm  and  1pm  x  20pm,  (c)  1pm  x  85  pm  and 
4pm  x  60pm,  and  (d)  1pm  x  100pm  and  1pm  x  130pm. 
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To  illustrate  different  trends  observed  in  Zn,  the  cross-coupling  impedances 
obtained  from  EPIC  simulations  for  four  different  contact  geometries  is  shown  in  Figure 
3.10. 

In  the  case  of  misaligned  contacts,  Zn  does  not  exhibit  a  consistent  exponential 
decrease  with  respect  to  GMD,  especially  for  long  contacts  (L  >  50pm)  at  close 
separations.  The  Zn  plot  with  Xee  =  2pm  portrays  near  field  effects  where  Zn  is  not  an 
exponential  function  of  the  GMD  for  long  contacts.  The  Zn  plot  when  Xee  is  very  large 
(50pm),  exhibits  far  field  effects  where  Zn  is  an  exponential  function  with  respect  to  the 
GMD. 

From  the  above  observations,  the  behavior  of  Zn  is  differentiated  based  on 
relative  positions  of  contacts  in  the  following  two  ways:  A  contact  within  the  proximity 
of  another  contact  and  a  contact  not  in  the  proximity  of  the  other  contact.  Figure  3.11 
shows  these  definitions. 


Proximity  regions  No-proximity  regions 


(a)  (b) 


Figure  3.11.  Relative  positions  of  Contact  2  with  respect  to  Contact  1  (a)  Contact  2  in  the 
proximity  region  of  Contact  1.  (b)  Contact  2  in  the  no-proximity  region  of  Contactl. 
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The  following  is  a  summary  based  on  an  analysis  of  the  trends  in  the  cross¬ 
coupling  impedances  for  the  entire  geometric  data  set: 

i.  The  cross-coupling  impedance  is  exponential  with  respect  to  the  GMD  for 
all  cases  when  the  second  contact  lies  out  of  the  proximity  region  of  the 
first  contact.  This  region  is  portrayed  in  Figure  3.1 1(b). 

ii.  The  deviation  from  the  exponential  behavior  is  prominent  for  long 
contacts  (L  >  50pm)  that  are  within  the  proximity  region  as  marked  in 
Figure  3.1 1(a). 

3.4.2. 1.  Z12  model  for  contact  cases  that  lie  outside  each  other’s  proximity  region 

Equation  (3.5)  developed  in  the  previous  section  is  applied  to  model  the  mutual 
impedance  between  contact  geometries  that  are  located  outside  each  contact’s  proximity 
region  (Figure  3.12). 


Positions  of 
Contact  2 


Figure  3.12.  Positions  for  Contact  2  with  respect  to  Contact  1  that  are  outside  the 
proximity  region  of  Contact  1 . 
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The  relative  errors  between  the  EPIC  simulation  results  and  the  GMD  model  are 
within  ±10%  for  the  same  range  of  contact  dimensions  as  quantified  in  Section  3.4.1 
(Figures  3.13  and  3.14). 


Y>m>  Yc>m) 

(C)  (d) 


Figure  3.13.  Comparison  of  the  GMD  model  with  EPIC  simulations  for  cases  placed  in 
the  no-proximity  region  withXee  =  4pm.  The  contact  dimensions  are  (a)  5pmx5pm  and 
5pmx5pm,  (b)  20pmxlpm  and  lpmx20pm,  (c)  lpmx85pm  and  4pmx60pm,  and  (d) 

lpmxlOOpm  and  lpmxl30pm. 
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Figure  3.14.  Comparison  of  the  GMD  model  with  EPIC  simulations  for  cases  placed  in 
the  no-proximity  region  withXee  =  50pm.  The  contact  dimensions  are  (a)  5pmx5pm  and 
5pmx5pm,  (b)  20pmxlpm  and  lpmx20pm,  (c)  lpmx85pm  and  4pmx60pm,  and  (d) 

lpmx  100pm  and  lpmx  130pm. 


3.4.2.2.  Z12  model  for  contact  cases  that  lie  inside  the  proximity  region 

For  the  cases  where  a  contact  is  in  the  proximity  region  of  another  contact,  Zu  is 
not  an  exponential  function  of  the  GMD.  To  model  the  behavior  of  Z12  in  these  regions, 
an  analysis  was  carried  out  with  different  contact  geometries  and  orientations.  It  was 
found  that  the  Z12  is  a  function  of  not  only  the  GMD  but  also  the  overlap-length  inside 
the  proximity  region.  The  idea  of  overlap-length  is  made  clear  in  Figure  3.15. 
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Figure  3.15.  Overlap-length  definitions. 


Two  examples  are  considered,  to  support  the  above  finding.  In  each  example,  one 
contact  is  placed  in  the  center  of  the  die,  and  the  other  contact  is  separated  by  increasing 
and  decreasing  Ycc  with  a  constant  Xee  (Figure  3.16).  In  the  first  example,  the  length  of  the 
2nd  contact  is  small  compared  to  the  length  of  the  1st  contact.  The  overlap-length  changes 
only  when  the  2nd  contact  is  near  the  corners  of  the  first  one  (Figure  3.17(c)).  The  cross¬ 
coupling  impedance  Z[2  behaves  in  the  same  manner  (Figure  3.17(a)).  In  the  second  case, 
the  lengths  of  the  contacts  are  comparable,  and  the  overlap-length  decreases  linearly  with 
the  maximum  being  at  Ycc=  0  (Figure  3.17(d)).  Zi2  follows  the  same  trend  as  shown  in 
Figure  3.17(b). 
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Example  1:  lpmx  110pm  Example  2:  lpm  x  100pm 

and  5pm  x27pm  and  lpm  x  1 00  pm 

Figure  3.16.  Relative  positions  for  Contact  2  with  respect  to  Contact  1  as  Ycc  is  varied. 

(Figures  not  drawn  to  scale) 
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Figure  3.17.  The  behavior  of  Zn  and  the  overlap-length  when  Ycc  is  increased,  (a)  Zu 
from  EPIC  simulations  for  Example  1.  (b)  Zn  from  EPIC  simulations  for  Example  2.  (c) 
Corresponding  overlap-length  behavior  for  Example  1 .  (d)  Corresponding  overlap-length 

behavior  for  Example  2. 

Clearly  Z12  depends  on  the  overlap-length  in  the  proximity  region.  The  point  of 
maximum  overlap  is  at  Ycc= 0.  The  impedance  at  Ycc= 0  is  computed  using  Equation  (3.5), 
since  the  exponential  behavior  of  Zi2  is  guaranteed  when  Xee  is  increased,  keeping  Ycc  =  0 
(Section  3.4.1).  In  order  to  model  the  behavior  of  Z12  with  proximity  effects,  Equation 
(3.5)  is  modified  to  include  a  normalized  overlap-length  term. 
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where  c?  is  a  function  of  the  geometry  and  Zl2(C-c)  is  the  cross-coupling  impedance  from 
Equation  (3.5).  It  is  to  be  noted  that  at  Ycc=  0,  the  overlap  is  at  its  maximum  value  and 
overlap  length  equals  maximum  overlap  length.  Hence  Equation  (3.6)  reduces  to 
Equation  (3.5). 

Equation  (3.6)  models  the  behavior  of  Zn  in  cases  where  the  second  contact  is  in 
the  proximity  region  of  the  first  contact.  The  performance  of  the  improved  model 
(Equation  (3.6))  is  compared  with  the  original  GMD  model  (Equation  (3.5))  and  with 
EPIC  simulations.  Results  for  four  contact  cases  with  different  Xee  values  are  presented  in 
Figure  3.18. 

The  plots  confirm  the  earlier  observation  that  the  original  GMD  model  suffers 
from  accuracy  issues  (relative  errors  >  15%)  at  close  separations.  On  the  other  hand,  the 
improved  model  has  relative  errors  that  are  less  than  10%.  The  relative  errors  in  the 
improved  model  increase  forXee  less  than  2pm. 

There  is  a  discontinuity  in  the  model  equations  in  the  region  of  transition  where 
the  improved  model  is  switched  to  the  GMD  model.  The  switching  of  the  improved 
model  to  the  GMD  model  is  required  to  describe  the  impedances  in  the  proximity  region 
and  no-proximity  regions,  respectively.  Standard  linear  or  exponential  smoothing 
techniques  were  not  efficient  in  making  the  transition  smooth  between  the  two  model 
equations  and  the  possibility  of  using  more  advanced  smoothing  functions  remains  a 
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direction  for  future  research.  The  discontinuity  causes  more  errors  compared  to  the  GMD 
model  for  a  few  cases  as  shown  in  Figure  3.18(a). 


Figure  3.18.  Comparison  of  the  improved  Zi2  model  equation  with  model  Equation  (3.5) 
and  EPIC  for  contact  dimensions  (a)  lpmx  1 10pm  and  5pmx73pm  (Xee=2[im),  (b) 
150pmx  150pm  and  150pmx  150pm  (Xee  =4pm),  (c)  3pmx80pm  and  4pmx  100pm  (Xee 
=4pm),  and  (d)  lpmxl  10pm  and  5pmx20pm  (. Xee  =10pm). 


The  next  step  is  to  model  a  as  a  function  of  the  geometry  and  identify  the 
proximity  regions  where  the  improved  model  equation  is  applicable.  The  approach  to 
modeling  c/  is  similar  to  that  of  ai  in  Section  3.4.1.  From  the  experiments  above,  cj  is 
extracted  for  various  contact  cases  at  a  particular  value  of  Xee.  The  values  are  shown  in 
Figure  3.19(a).  For  contacts  having  small  lengths,  the  proximity  effects  are  negligible. 
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Hence,  c/  should  be  such  that  Equation  (3.6)  reduces  to  Equation  (3.5)  for  contacts  with 
small  lengths.  A  logarithmic  function  is  used  to  model  c?  to  capture  this  behavior. 


c  i  =  log 


^  r  sum  _of  _ lengths  ^  ^ 


4 

V  v 


(3.7) 


2 


sum  _of  _  widths 

where  k4  and  k5  are  process  constants  with  values  k4=0.002  (pm)'0  75  and  k5=0.25.  For 
contacts  with  small  lengths,  cj  reduces  to  0. 
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Figure  3.19.  Accuracy  of  the  cj  model  for  all  contact  geometries,  (a)  Comparison  of  the 
cj  model  with  the  improved  model  extracted  values,  (b)  Relative  errors  for  the  c/  model 
compared  with  the  improved  model  extracted  values. 
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The  accuracy  of  this  model  is  determined  by  comparing  the  a  values  computed 
from  Equation  (3.7)  with  the  values  extracted  while  developing  the  improved  model. 
Figure  3.19(b)  shows  that  the  relative  errors  for  this  ci  model  are  within  12%. 

Next  the  region  of  proximity  is  considered.  The  first  observation  is  that,  the 
proximity  region  (Figure  3.11(a))  is  larger  for  longer  contacts.  For  example,  a  set  of 
square  contacts  with  dimension  100pm  may  experience  proximity  effects  at  a  separation 
of  50pm.  On  the  other  hand,  contacts  with  dimension  5pm  may  not  experience  proximity 
effects  at  the  same  separation.  The  ratio  of  the  sum  of  the  lengths  (L  i+Lj)  to  the  shortest 
distance  between  the  edges  of  the  contacts,  Xee,  is  used  to  quantify  these  regions.  The 
relative  error  of  the  GMD  model  was  observed  to  be  higher  than  15%,  due  to  proximity 
effects,  when  this  ratio  was  greater  than  10.  Consequently,  the  improved  model  equation 
is  used  whenever  the  ratio  exceeds  10. 

The  complete  cross-coupling  impedance  model  formulation  is  given  by  the  set  of 
equations: 


_  -k3-y]  gmd 

12  —  &\€  if  contacts  are  not  within  the  proximity  region 


(3.8) 


-k^gmdc_ 


f 


overlap_length 


maxmum_overlap_length 


if  contacts  are  within  the  proximity 


region  (3.9) 

The  gmdc.c  is  the  geometric  mean  distance  between  contacts  when  their  centers 
are  aligned.  The  proximity  effects  in  modeling  Zn  occur  when  the  ratio  of  the  sum  of 
contact  lengths  (Li+L2)  to  the  shortest  distance  between  the  edges  of  the  contacts,  Xee  is 
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greater  than  10.  The  final  comparison  of  the  entire  macromodel  with  the  EPIC  simulated 
values  is  illustrated  in  Figure  3.20. 


Figure  3.20.  Comparison  of  Zn  macromodel  with  EPIC  simulation  results.  The  contact 
dimensions  are  (a)  lpmxl  10pm  and  5pmxl08pm,  (b)  70pmx70pm  and  70pmx70pm, 
(c)  6pmx  180pm  and  5pmx  150pm,  and  (d)  30pmx6pm  and  6pmx90pm. 


The  complete  cross-coupling  impedance  model  requires  a  total  of  5  process 
constants.  The  relative  errors  between  the  model  and  EPIC  simulations  are  within  ±15% 


for  a  large  set  of  dimensions  and  separations. 
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3.5.  Self-impedance  (Zn)  model 

The  self-impedance  models  presented  in  [7]  and  [9]  are  functions  of  the  area  and 
perimeter  of  a  contact  and  do  not  account  for  any  interaction  from  the  neighboring 
contacts.  At  small  separations,  the  self-impedance  of  a  contact  is  affected  by  the 
surrounding  contacts  [16].  Figure  3.21  shows  the  variation  of  the  self-impedance  with 
separation  for  four  different  square  contacts. 


Figure  3.21.  Examples  to  show  that  Zn  deviates  from  constant  values  for  small 
separations.  The  contact  sizes  are  (a)  5pmx5pm  and  5pmx5pm,  (b)  lOpmxlOpm  and 
lOpmxlOpm,  (c)  20pm  x20pm  and  20pmx20pm,  and  (d)  50pmx50pm  and 

50pmx50pm. 


45 


It  is  clear  from  the  plots  that  Zn  decreases  with  a  decrease  in  separation  due  to  the 
presence  of  the  other  contact.  As  the  separation  increases,  this  effect  diminishes,  and  Zn 
approaches  a  constant  value  which  is  a  function  of  the  contact’s  geometry. 

As  a  first  step,  the  Zn  (or  Z22)  are  modeled  assuming  there  is  no  influence  from 
the  surrounding  contacts.  This  is  consistent  with  the  assumption  that  Z-parameters  are 
independent  of  other  contacts  (Section  2.2).  EPIC  simulations  were  performed  for  the  set 
of  contact  geometries  and  self-impedance  values  were  extracted.  The  resulting  self¬ 
impedance  values  are  modeled  using  the  contact  area  and  perimeter  information  as: 

v  1 

zn=  — - : - : - “Z -  (3.10) 

a j  {area)  +  a2  {perimeter)  3  +  a  4 

where  aL  a3,  a3  and  a 4  are  process  constants.  The  parameter  a3  is  introduced  to  increase 
the  accuracy  of  previous  models  [7],  [9].  Without  a3  the  model  gives  errors  greater  than 
20%  and  as  high  as  40%  when  compared  to  EPIC  simulations  results.  However,  the 
improved  model  matches  the  simulation  results  with  relative  errors  less  than  8%,  as 
depicted  in  Figure  3.22.  The  process  constants  for  this  model  are: 
on  =9.5x1 0\m)'2fi4  a2=3xlO'6(pm)'°'5Q'1  a3=0.5 

a4=4.6xl0'4f2'1 

Proximity  effects  due  to  the  neighboring  contacts  for  the  self-impedance  are  described  in 
Chapter  4. 
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Figure  3.22.  Zn  model  as  a  function  of  the  contact  geometry,  (a)  Comparison  of  Zn 
model  with  EPIC  values,  (b)  Relative  errors  for  the  Zn  model  compared  with  EPIC 

values. 


3.6.  Summary  of  the  macromodel  development 

This  chapter  was  dedicated  to  the  analysis  and  modeling  of  Z-parameters  for  two 
contacts  in  a  lightly  doped  substrate.  The  GMD  model  was  improved  using  a  simple 
overlap-length  function  to  eliminate  the  inaccuracies  caused  by  the  GMD  at  close 
separations.  This  is  crucial  as  practical  layouts  (both  analog  and  digital)  consist  of  closely 
spaced  long  contacts.  The  self-impedance  model  was  also  improved  using  an  extra 
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process  parameter.  The  self-impedance  model  does  not  account  for  the  proximity  of 
neighboring  contacts.  The  model  equations  (Equations  (3.8),  (3.9),  (3.10))  developed  in 
this  chapter  require  nine  process  constants  and  have  a  reasonable  accuracy  for  a  very 
large  set  of  layout  geometries  and  separations. 

The  geometric  data  set  for  the  model  parameter  extraction  has  been  optimized  in 
Appendix  E.  The  optimization  of  the  data  set  is  required  for  efficient  extraction  of  model 
parameters  when  the  macromodel  equations  are  needed  for  different  processes. 
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4.  PANELING  APPROACH  TO  OVERCOME  MACROMODEL  LIMITATIONS 


In  the  previous  chapter,  the  lightly  doped  substrate  was  characterized  based  on 
two-port  Z-parameters.  The  self-impedance  was  assumed  to  be  independent  of  the 
surrounding  contacts.  However,  this  assumption  is  not  valid,  and  the  self-impedance  is 
affected  by  surrounding  contacts  placed  at  close  separations.  This  chapter  examines  a 
paneling  approach  for  incorporating  these  proximity  effects  in  the  self-impedance  model. 
In  addition,  the  paneling  approach  is  applied  with  the  GMD  model  for  modeling  mutual 
impedances.  This  is  to  check  if  the  GMD  limitations  are  overcome  without  using  the 
improved  model.  Once  proximity  effects  are  accounted  for  in  the  two-port  model,  the 
macromodel  with  paneling  can  be  easily  extended  for  multiple  contacts. 

4.1.  The  need  for  paneling 

By  definition,  the  Z  or  open-circuit  parameters  are  independent  of  the  neighboring 
contacts.  However,  proximity  or  presence  of  near-by  contacts  affects  the  self-impedance 
values.  Figure  4.1  illustrates  the  variation  in  the  self  impedance  for  closely  placed 
contacts.  The  plots  show  normalized  self-impedance  values  for  different  contact  sizes  for 
a  range  of  separations.  It  is  observed  that  small  contacts  exhibit  proximity  effects  at 
smaller  separations  and  large  contacts  for  much  larger  separations,  i.e.,  greater  than 
100pm  for  two  lOOpmx  100pm  contacts.  Therefore,  the  size  of  the  contacts  as  well  as  the 
separations  have  to  be  considered  for  proximity  effects.  For  this  reason,  the  previous  Z- 
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parameter  based  model  is  not  adequate  and  the  influence  of  the  neighboring  contacts  on 
Zn  (or  Z22)  has  to  be  incorporated. 

It  is  tedious  to  include  the  influence  of  neighboring  contacts  while  modeling  the 
self-impedances  for  an  N-contact  network.  The  impedance  of  an  N-port  is  represented  by 
an  NxN  impedance  matrix  Z.  Due  to  contact  proximity  each  diagonal  element  (Zn  or  self¬ 
impedance)  of  Z  will  depend  on  the  other  N-l  contacts.  This  complicates  the 
dependencies  of  Zn.  A  paneling  approach  has  been  suggested  in  [17]  to  model  such 
effects.  This  idea  was  based  on  the  work  presented  in  [18].  The  paneling  approach  is  also 
used  in  this  thesis.  A  ±5-10%  variation  in  the  self-impedance  values  is  assumed  to  be 
acceptable  for  small  contacts. 


Figure  4.1  Influence  of  proximity  effects  on  self-impedance  values. 
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4.2.  The  paneling  approach 

The  paneling  approach  involves  dividing  the  contacts  into  panels.  The  panels  are 
assumed  to  be  equipotential  and  have  a  unifonn  current  density.  The  paneling  algorithm 
has  been  developed  for  a  general  case  where  the  contacts  have  been  divided  along  both 
the  length  and  width  as  shown  in  Figure  4.2. 


^  w 

Separation 

Contact  1  divided  into  a  panels  Contact  2  divided  into  b  panels 

Figure  4.2.  Panel  division  for  2  contacts. 

The  paneling  algorithm  steps  are  as  follows: 

Step  1 :  Panel  formation 

•  Divide  the  contacts  into  a  and  b  panels,  respectively. 

•  Obtain  the  sizes  and  separations  of  panels  based  on  geometric  information 
to  calculate  the  panel  level  Z-matrix  (Zp). 

Step  2:  Calculating  the  panel  level  Z-matrix  (Zp) 

The  calculation  of  the  panel  level  ( a+b)x(a+b )  Z-matrix  is  shown  in  Figure  4.3. 
Based  on  the  relative  positions  of  panels,  the  panel  level  Z-parameters  are  calculated 
using  the  macromodel  equations  (Equation  (3.8),  Equation  (3.9),  and  Equation  (3.10)). 
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Figure  4.3.  Flowchart  for  Zp  calculation. 


Step  3:  Calculation  of  the  2-port  Z-matrix 


•  Yp  =  inv(Zp);  The  structure  of  this  admittance  matrix  is  as  follows: 


Y„ 


•  Calculating  the  2-port  R-matrix 

a.  The  self  resistance  is  calculated  as  below 
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b.  The  mutual  resistance  is  calculated  as: 
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The  final  R-matrix  is: 
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•  Calculate  2-port  Y-matrix  from  the  above  2-port  R-matrix 

•  Calculate  2-port  Z-matrix  by  inverting  the  2-port  Y-matrix 

Step-2  in  the  above  approach  requires  that  the  macromodel  be  applied  to  adjacent  panels. 
Before  proceeding  to  validate  our  paneling  approach,  an  explanation  of  how  the 
macromodel  can  be  applied  to  adjacent  panels  is  provided. 

4.2.1,  The  impedance  model  for  adjacent  panels  (or  contacts) 

Since  the  two  panels  are  adjacent  to  each  other,  they  can  be  considered  to  be 
equipotential.  Two  adjacent  contacts  are  represented  in  Figure  4.4(a)  and  the  merged 
single  contact  is  shown  in  Figure  4.4(b).  The  equivalent  n -resistive  network  is 
represented  in  Figure  4.4(c). 

adjacent  contacts  merged  contact 


(c) 

Figure  4.4.  (a)  Adjacent  contacts,  (b)  Merged  contact,  (c)  Resistive  network  for 

equipotential  adjacent  contacts. 
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The  cross-coupling  resistance,  Rij,  for  the  adjacent  contact  network  is  zero. 
Therefore,  the  two-port  network  reduces  to  a  single  port  network: 


z  RaRjj  _  RaRjj 

"  Ru  +  Rjj  +  R{j  Ru+Rjj 


merged  _  contacts 


(4.1) 


and 


Zii  ~  Rii  II  (Rjj  +  Rij)  ~  Rii  II  Rjj 


merged  _contacts 


(4.2) 


Zjj=Zii=Z=self-impedance  of  the  overall  merged  contact 


EPIC  was  used  to  perfonn  simulations  to  detennine  the  validity  of  the  above 
equations,  i.e.,  can  two  adjacent  contacts  (or  panels)  be  considered  as  a  merged  contact 
(Figure  4.4).  Many  different  geometric  cases  have  been  simulated,  each  with  contacts 
having  different  dimensions.  The  separation  between  the  adjacent  contacts  is  varied 
between  Opm  and  1pm.  The  results  for  two  test  cases  are  presented  in  Figure  4.5.  It  can 
be  observed  from  the  plots  that  the  circuit  parameters  for  adjacent  contacts  (Z12,  Zn)  do 
not  match  the  self-impedance  of  the  merged  contact,  Zmerged  (at  separation=0pm),  which 
is  not  the  expected  behavior  for  adjacent  contacts.  EPIC  does  not  consider  two  adjacent 
contacts  as  equivalent  to  a  single  merged  contact  contrary  to  Equations  (4.1)  and  (4.2). 
This  is  clearly  a  problem  with  EPIC.  For  this  reason,  the  adjacent  panels  or  contacts  are 
not  modeled  as  a  merged  single  contact.  Instead,  a  2-port  network  is  extracted  based  on 
adjacent  contacts  being  treated  as  very  closely  spaced  contacts. 
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Figure  4.5.  Comparison  of  Z-parameters  from  EPIC  for  merged  single  contacts  with 
adjacent  two  contacts  having  contact  dimensions  (a)  lpmxlpm  and  1  umd  urn,  and  (b) 

5pmx5pm  and  5pmx5pm. 


Hence,  the  cross-coupling  impedance  for  adjacent  contacts  can  be  approximated 
by  Zi2A-  Z]2a  denotes  the  cross-coupling  impedance  for  2  contacts  placed  at  a  very  small 
separation  (~  0.5pm).  The  self  impedances  are  given  by  Equation  (3.10).  The 
macromodel  using  small  separations  to  model  adjacent  contacts,  is  shown  to  agree  with 
relative  errors  less  than  ±15%  with  EPIC  values.  The  results  have  been  tabulated  in  Table 


4.1  for  a  number  of  cases. 
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Table  4.1.  Comparison  of  Zn  values  from  EPIC  and  the  macromodel  for  very 


closely  spaced  panels/contacts  (~  0.5pm). 


Contact  Geometry 

EPIC  values  (f2) 

Macromodel  (f2) 

Error  % 

6 pm x 6 pm  and  6pm><6pm 

451.19 

385.73 

14.5% 

10pmx3pmand  3pmxl6pm 

407.63 

366.44 

10.0% 

lOpmx  10pm  and  10pm  x  10pm 

380.78 

342.22 

10.0% 

125pmxlpmand  125pmxlpm 

115.13 

108.76 

5.5% 

lpmx30pmand  2pmx2pm 

403.85 

394.98 

2.2% 

40pmx40pm  and  40pmx40pm 

204.09 

202.8 

0.063% 

4.2.2.  Paneling  technique 

While  developing  the  paneling  approach,  it  is  important  to  establish  the  direction 
in  which  paneling  is  performed  along  the  width  of  a  contact,  length  of  a  contact,  or  both. 
Experiments  to  determine  which  paneling  technique  provides  a  better  estimate  of  self¬ 
impedances  have  been  carried  out.  Consider  a  pair  of  square  contacts  placed  relative  to 
each  other,  such  that  neighboring  contacts  affect  the  self-impedances  (Section  4.1).  The 
paneling  algorithm  was  applied  to  calculate  the  Z-parameters.  EPIC  simulations  were 
used  to  calculate  the  Zp  matrix  (Section  4.2),  to  accurately  predict  the  self-impedances. 
The  self-impedances  from  the  paneling  approach  were  compared  to  the  self-impedances 
extracted  directly  from  EPIC. 

Paneling  was  carried  out  in  two  different  ways.  The  contacts  were  divided  into 
panels  along  the  length  (vertical  paneling)  and  along  the  width  (horizontal  paneling). 
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Figure  4.6(b)  and  Figure  4.7(b)  show  that  for  the  given  pair  of  contacts,  vertical  paneling 
was  more  effective  compared  to  horizontal  paneling. 


(a)  (b) 

Figure  4.6.  Effect  of  vertical  paneling,  (a)  Paneling  technique,  (b)  Comparison  of  Zn 
obtained  from  varying  panel  sizes  with  reference  Zu,  for  a  pair  of  20pmx20pm  contacts. 


(a)  (b) 

Figure  4.7.  Effect  of  horizontal  paneling,  (a)  Paneling  technique,  (b)  Comparison  of  Zn 
obtained  from  varying  panel  sizes  with  reference  Zn,  for  a  pair  of  20umx20um  contacts. 
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Next  consider  two  long  skinny  contacts.  The  paneling  approach  is  applied  to  this 
example.  EPIC  simulations  were  carried  out  to  extract  the  self-impedances.  The  results 
for  horizontal  and  vertical  paneling  are  shown  in  Figure  4.8.  It  was  found  that  horizontal 
paneling  was  more  effective  compared  to  vertical  paneling. 


Figure  4.8.  Horizontal  and  vertical  paneling  for  long  skinny  contacts  (3pmx80pm  and 

4pmx  100pm). 


Since  the  accuracy  of  the  paneling  technique  is  geometry  dependent,  simultaneous 
horizontal  and  vertical  paneling  was  adopted.  This  makes  the  paneling  algorithm  more 
generic  and  applicable  to  a  wide  variety  of  contact  geometries  (Figure  4.2). 
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4.3.  Macromodel  with  paneling  for  two-contact  cases 

The  macromodel  along  with  the  paneling  approach  described  in  the  previous 
section  is  applied  to  two-contact  examples  in  this  section.  This  is  done  to  check  if  the 
proximity  effects  on  the  self-impedances  are  properly  accounted  for.  Also,  the  paneling 
approach  is  applied  to  model  mutual  impedances  using  the  GMD  model  to  examine  if  the 
paneling  approach  can  overcome  the  GMD  limitations. 

4.3.1.  Self-impedance  model  with  paneling 

A  simple  example  of  a  pair  of  square  contacts  is  considered  and  the  paneling 
approach  is  verified  for  different  panel  sizes.  The  modeled  self  impedances  have  been 
compared  with  the  extracted  impedance  values  from  EPIC.  Figure  4.9  shows  that  for 
panel  sizes  less  than  5  pm,  Zn  modeled  by  the  paneling  approach  does  not  tend  towards 
the  actual  Zn  values  extracted  from  EPIC. 


Figure  4.9.  Macromodel  with  paneling  applied  to  a  set  of  25pmx25pm  contacts. 
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If  the  contacts  are  divided  into  more  number  of  panels,  the  resulting  Zn  should  be 
more  accurate.  However,  deviations  from  this  expectation  can  be  attributed  to  the 
inaccuracies  in  the  macromodel,  (for  small  dimensions  at  close  separations)  as  stated  in 
Section  3.4.  Since  the  proximity  effects  are  insignificant  for  small  contacts,  the  paneling 
approach  is  applied  with  a  minimum  panel  size  of  5pm.  The  paneling  approach  was  also 
applied  to  larger  contacts,  and  the  proximity  effects  were  properly  modeled.  The  results 
for  four  different  contact  geometries  are  shown  in  Figure  4.10. 


Separation  (^m)  Separation  (urn) 

(c)  (d) 

Figure  4. 10.  Macromodel  Zn  with  paneling  compared  with  EPIC  Zn  for  contact 
dimensions  (a)  50pmx50pm  and  50pmx50pm,  (b)  30pmx6pm  and  6pmx90pm,  (c) 
125pmxlpm  and  125pmxlpm,  and  (d)  25pmx20pm  and  16pmx7pm. 
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The  relative  errors  were  observed  to  be  less  than  5%  when  compared  to  EPIC 
simulations.  Clearly  the  macromodel  applied  without  paneling  gives  errors  greater  than 
15%  when  compared  to  EPIC  simulations. 

4.3.2  Mutual-impedance  model  with  paneling 

In  the  previous  chapter,  the  GMD  limitations  for  the  mutual  impedance  model 
were  overcome  using  an  improved  model  (Equation  (3.9)).  In  this  chapter,  the  paneling 
approach  is  used  in  conjunction  with  the  GMD  model  (Equation  (3.8))  to  check  if  the 
GMD  limitations  are  resolved.  The  results  for  two  examples  are  shown  in  Figure  4.1 1. 


Figure  4.1 1.  Macromodel  Z12  with  paneling  compared  with  improved  model  and  EPIC 
Z12  for  contact  dimensions  (a)  lpmxl  10pm  and  5pm><73pm  placed  atXee=2pm,  and  (b) 
150pmx  150pm  and  150pmx  150pm  placed  atXee=4pm. 

Other  geometric  cases  and  proximity  regions  that  required  the  application  of  the 
improved  model  were  also  tested  with  the  GMD  model  in  conjunction  with  paneling.  It 
was  found  that  the  GMD  model  with  paneling  performed  better  than  the  improved  model 
at  the  expense  of  computation  time.  For  example,  to  acquire  one  value  of  Zi2  for  the  pair 
of  150pmxl50pm  contacts,  the  macromodel  with  paneling  (panel  size=10pm)  took  20 
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times  more  computation  time  than  the  macromodel  without  paneling.  However,  because 
of  the  paneling  approach,  relative  errors  for  all  cases  were  of  the  order  of  10%  when 
compared  to  EPIC  simulations.  There  were  no  issues  of  discontinuity  with  this  approach, 
when  paneling  is  applied  irrespective  of  the  proximity  region. 

4.4.  Extension  of  the  model  to  multi-contact  cases 

As  a  first  step  in  extending  the  model  to  multi-contact  examples,  a  three-contact 
case  is  considered  (Figure  4.12).  Because  of  the  presence  of  a  third  contact,  the  Z- 
parameters  are  affected  by  the  neighboring  contacts.  The  three  contact  substrate  network 
is  modeled  using  the  macromodel  in  conjunction  with  paneling  to  account  for  the 
proximity  of  the  third  contact. 

Examples  of  three  contact  Z-parameters  are  shown  in  Figure  4.13.  Both  the  Zn 
and  Zjj  are  modeled  with  reasonable  accuracy  using  the  macromodel  with  paneling.  The 
errors  are  within  10%  for  both  the  self  and  mutual  impedances  with  this  approach.  The 
macromodel  without  paneling  for  the  three  contact  example  has  errors  as  high  as  60%. 


50  pm 


Ycc  Ycc 


50  pm  5  pm  50  pm 


Figure  4.12.  Example  of  3  contacts  separated  by  Ycc distances. 
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+  EPIC,  - Macromodel  with  paneling,  +  Macromodel  without  paneling 


Figure  4.13.  Comparison  of  the  macromodel  with  paneling  with  EPIC  simulations  for  a 

three  contact  example. 
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Another  example  of  a  multi-contact  case  is  considered  to  evaluate  the 
performance  of  the  macromodel  applied  with  paneling.  Five  square  contacts  with 
dimensions  50pm  were  placed  on  a  Immxlmm  die  as  shown  in  Figure  4.14.  The  Xee 
separation  was  varied  and  the  Z-parameters  were  extracted  using  EPIC. 

600 

550 

W 

I  500 
>- 

450 

400 

400  450  500  550  600 

X  axis 

Figure  4.14.  Layout  of  five  square  contacts  with  dimension  50qm  in  a  Immxlmm  die. 


5  contacts 


T - 1 - ' - 1 


J _ ! _ L 


(a) 


(b) 


Figure  4.15.  Comparison  of  the  macromodel  with  paneling  and  macromodel  without 
paneling  with  EPIC  simulations  for  a  five  contact  example,  (a)  Z12.  (b)  Z\\. 
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Figure  4.16  presents  the  results  for  Zn  and  Zn.  The  comparison  between 
macromodel  (with  and  without  paneling)  and  EPIC  simulations  are  shown.  It  can  be 
observed  that  the  paneling  approach  is  required  to  model  the  multi-contact  Z-parameters 
correctly. 

Based  on  the  above  analysis,  it  is  clear  that  the  paneling  approach  is  important  for 
extending  the  two-contact  model  to  multiple  contacts.  In  Chapter  5,  a  complex  multi¬ 
contact  example  is  validated  using  this  extension.  An  outline  of  the  advantages  and 
pitfalls  of  the  mutual-impedance  modeling  approaches  is  given  in  Table  4.2.  The  extra 
computation  time  for  the  paneling  approach  depends  on  the  contact  geometry  and  the 
panel  sizes.  A  summary  of  the  performance  of  the  entire  macromodel  for  two-contact  and 
multi-contact  cases  is  listed  in  Table  4.3. 


Table  4.2.  Comparison  of  improved  Zn  model  with  the  paneling  approach. 


Aspect 

Improved  Z  n  model 

GMD  model  with  paneling 

Relative  errors 
compared  to  EPIC 

within  ±10-15% 

within  ±10% 

Computation  time 

(100  values) 

<  2sec  (both  cases) 

•  4  min  (50pmx50pm  and  50pmx50pm 
for  a  minimum  panel  size=  l  0pm) 

•  20  sec  (lpmx  110pm  and  5pmx73pm 
for  a  minimum  panel  size=5pm  ) 

Discontinuity 

Yes 

No 
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Table  4.3.  Summary  of  the  model  performances  for  two-contact  and  multi-contact 

cases. 


Z-  Models 

2  Contacts 

Multiple 

Contacts 

w/o  paneling 

with  paneling 

w/o  paneling 

with  paneling 

Zi2:  GMD 
Model 

Eq.  (3.8) 

Errors  within 
±10%  when 
contacts  out  of 
proximity 
region. 

Errors  >30% 
for  contacts  in 
proximity 
region. 

Errors  within 
±10%  for  all 

cases 

>20%  errors 

Errors  within 
±10%  for  all 

cases 

Zj2:  Improved 
Model 

Eq.  (3.9) 

Errors  within 
±15%  for 
contacts  in 
proximity 
region. 

Not  used* 

>20%  errors 

Not  used* 

Zn:  Geometric 
Model 

Eq.  (3.10) 

Errors  within 
±8%  for 
contacts 
having  no 
proximity 
effects. 

>20%  errors 
for  contacts 
having 
proximity 
effects. 

Errors  within 
±5%  for  all 

cases 

>20%  errors 

Errors  within 
±5%  for  all 

cases 

*  The  paneling  approach  uses  a  uniform  panel  size  of  5pm.  Hence,  if  the  macromodel  is  used  with  paneling, 
the  improved  model  is  not  required. 
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5.  VALIDATION  OF  THE  MACROMODEL  FOR  LIGHTLY  DOPED 

SUBSTRATES 

The  macromodel  developed  in  Chapter  3  is  further  validated  for  a  different  set  of 
geometries  and  variations  of  the  contact  separation  along  the  diagonal.  Comparisons  with 
existing  measured  Z-parameter  values  obtained  for  a  0.25pm  CMOS  process  are  also 
made.  The  application  of  the  model  (in  conjunction  with  paneling)  for  simulating  the 
substrate  noise  coupled  to  a  2.4GHz  RF  LNA  from  a  stepped  buffer  is  also  presented. 

5.1.  Validation  for  a  different  data  set 

To  further  validate  the  macromodel,  it  is  essential  to  show  that  it  works  for  a  data 
set  that  that  was  not  considered  during  the  model  development  phase. 

Figure  5.1  presents  sixteen  test  cases  for  validating  the  self-impedance  model. 
The  relative  errors  between  the  macromodel  and  EPIC  simulations  are  within  ±5%. 
Figure  5.2  compares  the  mutual  impedance  extracted  from  EPIC  and  the  values  predicted 
by  the  model  for  four  different  contact  geometries.  The  model  is  in  very  good  agreement 
with  the  EPIC  extracted  Z-parameters  with  relative  errors  within  10%. 


Relative  Error  % 


68 


(a) 


(b) 


Figure  5.1.  Validation  of  the  self-impedance  macromodel  for  different  geometries,  (a) 
Comparison  of  the  Z\\ model  with  EPIC  values,  (b)  Relative  errors  for  the  Zn  model 

compared  with  EPIC  values. 
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Figure  5.2.  Validation  of  the  cross-coupling  impedance  macromodel  using  different 
geometries  for  Xee  and  Ycc  with  contacts  dimensions  (a)  7.5pmx7.5pm  and 
7.5pmx7.5pm,  (b)  lOpmxlOpmand  45  pmx45pm,  (c)  lpmx45pmand  5pmx50pm,  and 

(d)  lpmx93pmand  lpmx93pm. 

5.2.  Validation  of  the  model  for  diagonal  separation 

Next,  we  validate  the  model  by  varying  the  diagonal  separation  between  the 
contacts  by  simultaneously  changing  Xee  and  Ycc  as  shown  in  Figure  5.3.  In  Figure  5.4  the 
Z-parameters  extracted  from  EPIC  and  the  values  predicted  by  the  model  have  been 
compared.  The  Zu  model  is  found  to  be  in  good  agreement  with  the  EPIC  values. 
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Figure  5.3.  Contacts  separated  along  the  diagonal. 


(c)  (d) 


Figure  5.4.  Comparison  of  model  with  EPIC  simulations  for  diagonal  separation  with 
contact  dimensions  (a)  18pmx5pm  and  2pmxl5pm,  (b)  30pmx30pm  and  30pm  x30pm, 
(c)  1pm  xlpm  and  5pmx73pm,  and  (d)  lpm  xl  10pm  and  5pmx  100pm. 
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5.3.  Validation  using  measurement  results 

The  next  set  of  validations  for  the  macromodel,  use  resistance  measurements  for  a 
TSMC  0.25pm  lightly  doped  substrate.  Figure  5.5  shows  a  cross  section  of  the 
approximate  doping  profile  for  the  substrate.  While  applying  the  process  constants 
(obtained  for  the  typical  0.25pm  process)  from  Chapter  3,  to  model  the  self  and  mutual 
impedances  for  this  substrate,  the  errors  were  observed  to  be  as  high  as  40%. 


p+  channel-stop  implant 


p-type  substrate 


Figure  5.5.  Layered  doping  profile  for  a  lightly  doped  substrate. 


Hence,  measurements  of  9  test  structures  (chosen  randomly)  were  used  for 
estimating  the  mutual-impedance  process  parameters  and  2  were  used  for  validating  the 
model.  The  test  structures  haveXee  separations  varying  from  1pm  to  100pm.  So  Equation 
(3.5)  alone  is  used.  The  extracted  process  parameters  for  the  mutual-impedance  are: 
ki=2.22xlO'14Q  (pm)'1  k2=  611.15ft  k3=0.158  (pm)'05 

The  validation  results  are  shown  in  Figure  5.6.  The  relative  errors  in  cross-coupling 
impedances  are  about  15-20%.  The  errors  are  slightly  higher  than  expected  since  the 
model  parameters  were  extracted  from  a  limited  number  of  test  structures. 
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Figure  5.6.  Comparison  of  the  mutual  impedance  macromodel  with  measurement  values 
for  contact  dimensions  (a)  0.66pmx0.62pm  and  0.66pmxl.l8pm,  and  (b) 
0.66pmx0.62pm  and  0.66pmxl  1.98pm. 


For  the  self-impedance  model,  all  the  test  structures  were  used  to  extract  the 
process  constants.  The  results  are  shown  in  Figure  5.7.  The  errors  relative  to 
measurements  are  within  ±6%  and  the  process  parameters  were  calculated  as: 
ai  =9.5xlO'8(pm)'2fT1 
a4=5.56xl0'5f2'1 


a2=3xlO'6(pmy0-5f2-1 


a3=0.377 
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Figure  5.7.  Validation  of  the  self-impedance  macromodel  with  measured  values. 

(a)  Comparison  of  macromodel  with  measurements,  (b)  Errors  relative  to  the 

measurements. 

5.4.  Multi-contact  Problem 

The  model  has  also  been  applied  to  a  multi-contact  example.  A  complex  post 
layout  example  of  an  RF  LNA  is  used.  An  analysis  of  the  substrate  noise  for  the  LNA  has 
been  presented  in  [19].  A  stepped  buffer  was  used  as  the  noise  injector.  The  block 
diagram  for  the  stepped  buffer  and  a  circuit  diagram  for  the  LNA  are  provided  in  Figure 
5.8.  The  operating  frequency  was  set  at  2.4GHz  and  hence  the  substrate  network  was 
modeled  as  a  resistive  network.  Figure  5.9  shows  a  simplified  layout  of  the  circuits.  The 
network  consists  of  5 1  contacts. 

A  simple  model  depicting  the  substrate  network  is  shown  in  Figure  5.10.  Noise  is 
injected  into  the  substrate  because  of  the  switching  activity  of  the  inverters  in  the  stepped 
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buffer.  The  injected  noise  propagates  via  the  substrate  and  is  picked  up  by  the  sensitive 
LNA. 


EPIC  was  used  to  extract  the  equivalent  substrate  network  in  [18].  In  this  work, 
the  macromodel  in  conjunction  with  the  paneling  approach  is  used  to  extract  the 
equivalent  substrate  network. 

Stepped  buffer 


Vdd 


Input 


Figure  5.8.  Block  diagram  for  the  stepped  buffer  and  circuit  schematic  of  the  LNA. 
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This  section  evaluates  the  macromodel  performance  in  terms  of  computation  time 
and  accuracy  in  comparison  with  EPIC  for  the  circuit  example.  In  Appendix  F,  a  more 
complicated  substrate  network  for  the  same  circuit  example  (84  contacts)  is  considered 
for  comparison  of  the  macromodel  with  measurement  results. 
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Figure  5.9.  Fayout  of  the  FNA  and  stepped  buffer  with  5 1  contacts. 
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Figure  5.10.  Cross  section  of  the  substrate  network  for  LNA  and  stepped  buffer  in  a 

lightly  doped  substrate. 

The  frequency  spectrum  of  the  LNA  output  is  analyzed  to  evaluate  the  substrate 
noise.  The  important  resistance  values  have  been  compared  in  Table  5.1.  Though  some  of 
the  errors  are  of  the  order  of  20%  they  do  not  affect  the  output  noise.  The  model  applied 
with  paneling  is  in  good  agreement  with  EPIC  simulations  as  shown  in  Figure  5.10.  The 
model  applied  without  paneling  does  not  predict  the  output  noise  spectrum  as  close  as  the 
model  applied  with  paneling.  The  computation  time  for  the  model  with  paneling  is  only  9 
minutes  and  30  seconds,  when  compared  to  EPIC  that  took  5  hours  and  38  minutes  on  a 
Redhat  Linux  (kernel  version  2. 6. 9-34. EL)  machine  with  1GB  RAM  and  a  3.2GHz  CPU. 
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Table  5.1  Comparison  of  resistance  values  from  EPIC  and  those  generated  from  the 


macromodel. 


Resistance  values 

EPIC  (kO) 

Macromodel  (kO) 

Error  % 

Ru  (LNA  input  transistor) 

36.626 

33.571 

8.34% 

Rjj  (LNA  output  transistor) 

35.08 

39.43 

-12.4% 

Rkk  (Buffer  active  region) 

1.372 

1.5482 

12% 

Rko  (Buffer  active  region  to 

ground) 

0.025 

0.030 

20% 

Ru  (Buffer  ground  taps) 

0.561 

0.515 

8.3% 

Rnn  (Interconnect) 

0.349 

0.339 

2.8% 

2.2 


2.25 


2.3 


2.35  2.4  2.45 

Frequency  (GHz) 


2.5 


2.55 


xIO 


2.6 

.9 


Figure  5.1 1.  Comparison  of  spectrum  of  the  LNA  output  using  EPIC  and  the  macromodel 

for  the  substrate  network. 
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6.  CONCLUSIONS 


Substrate  noise  coupling  is  a  critical  issue  when  integrating  analog,  RF,  and 
digital  circuits  on  a  common  substrate.  Lightly  doped  substrates  are  preferred  for  analog 
and  RF  applications  as  they  provide  better  noise  isolation.  Efficient  macromodels  are 
needed  for  estimating  the  noise  coupling  in  large  systems-on-a-chip. 

This  work  has  presented  a  Z-parameter  based  macromodel  that  characterizes  the 
lightly  doped  substrate  as  a  resistive  network.  The  model  is  scalable  with  contact 
geometry  and  separations.  The  two-contact  model  requires  a  total  of  nine  process 
constants.  The  relative  errors  between  the  macromodel  and  EPIC  simulation  values  are 
about  10-15%  for  a  large  set  of  contact  dimensions  and  separations. 

The  cross-coupling  impedance  has  been  modeled  using  an  improved  geometric 
mean  distance  model.  The  modification  alleviates  the  inaccuracies  in  the  original  GMD 
model  for  close  separations.  This  is  important  as  such  cases  are  often  encountered  in 
practical  layouts.  The  self-impedance  model  is  enhanced  with  a  paneling  approach  that 
accounts  for  the  proximity  of  neighboring  contacts.  The  paneling  approach  can  also  be 
used  for  resolving  the  GMD  limitations  in  the  mutual-impedance  model  at  the  expense  of 
computation  time.  The  entire  model  is  in  good  agreement  with  simulation  results. 

The  backplane  of  the  lightly  doped  substrate  has  been  approximated  as  a  single 
node  in  this  work.  This  assumption  has  been  validated  for  thick  substrates  exceeding 
200pm.  The  relationship  between  the  die  area  and  the  substrate  thickness  for  this 
assumption  to  be  valid  is  not  explored.  This  aspect  needs  to  be  examined  further. 
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The  thesis  focuses  on  modeling  the  resistive  substrate.  With  an  increase  in  the 
frequency  of  operation,  the  substrate  network  becomes  capacitive  [20].  This  frequency 
dependence  needs  to  be  modeled,  for  extending  the  application  of  the  macromodel  to 
higher  operating  frequencies. 
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APPENDIX  A 

EPIC  (Extraction  of  Parasitics  in  Integrated  Circuits)  is  a  Green’s  function-based 
solver.  The  tool  characterizes  the  silicon  substrate  using  a  multilayered  structure.  The 
resistivity  and  thickness  of  each  of  these  layers  is  obtained  from  the  substrate  doping 
profile.  Solving  the  Green’s  function  yields  the  potential  at  a  point  on  the  substrate  due  to 
an  electric  charge  at  another  point.  For  more  information  on  the  solver  options,  refer  to 
[5]. 

The  version  of  the  tool  used  for  this  work  characterizes  the  substrate  as  a  resistive 
network  assuming  low  frequency  operation.  Also,  the  tool  assumes  contacts  on  the 
substrate  to  be  2-D,  i.e.,  there  is  no  depth  specification  required  for  describing  the 
contacts.  In  addition,  adjacent  contacts  are  not  automatically  considered  physically 
merged.  The  input  data  has  to  specify  adjacent  contacts  to  be  equipotential. 

The  discrete  cosine  transforms  (DCT)  are  used  by  EPIC  to  perfonn  a  fast 
computation  of  the  Z-matrix.  By  default,  the  backplane  of  a  substrate  is  grounded  by 
adding  a  low  resistive  layer  to  the  substrate’s  backside.  To  keep  the  backplane  floating, 
the  input  file  has  to  be  modified  by  adding  a  high  resistive  layer  to  the  backside  of  the 


substrate. 
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APPENDIX  B 

To  show  that  the  backplane  of  a  lightly  doped  silicon  substrate  can  be  considered 
an  equipotential  surface  when  the  substrate  thickness  is  large,  a  multi-contact  example  is 
considered.  The  contacts  have  large  dimensions  and  are  placed  in  a  1  nunx  1mm  die.  The 
layout  is  shown  in  Figure  B-l.  Table  B-l  compares  some  of  the  effective  coupling 
resistance  values  with  worst  case  errors  when  the  backplane  is  assumed  equipotential  for 
a  substrate  thickness  of  250  pm. 
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Figure  B-l  Layout  of  five  big  contacts  in  a  Immxlmm  die. 
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Table  B-l.  Comparison  o f  effective  cross-coupling  resistances  from  a  floating  backplane 
setup  and  those  derived  from  a  grounded  backplane  for  a  multi-contact  system  in  a 
_ _ 250pm  thick  substrate. _ _ 


Rij 

From  floating 

Extracted  from  grounded 

Relative  Error 

backplane(f2) 

backplane(fl) 

Rl2 

212.65 

206 

3% 

Rl4 

290.45 

271.6 

6.4% 

R25 

334 

280 

16% 

R45 

283.35 

260.09 

8% 

The  above  multi-contact  example  is  now  tested  using  a  resistive  ground 
connection.  Instead  of  directly  grounding  the  substrate  backplane,  a  resistive  layer  is 
added  to  simulate  an  improper  ground  connection  (Figure  B-2). 


p+  channel-stop  implant 


p-type  substrate 


Layer  to  provide  small 
resistance  from  the  backplane 
to  ground  connection 


Figure  B-2.  Cross  section  of  the  doping  profde  for  a  typical  0.25pm  process  lightly  doped 

substrate  with  a  resistive  ground  connection. 


It  was  observed  that  the  equipotential  backplane  assumption  holds  good  even  with 
a  resistive  ground  connection.  In  fact  there  were  no  changes  in  the  effective  cross- 
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coupling  resistance  values  derived  from  the  two  situations:  grounding  the  substrate 
backplane  directly  and  grounding  the  substrate  backplane  through  a  small  resistance. 

For  a  substrate  thickness  of  100pm,  the  multi-contact  example  is  used  for  testing 
the  equipotential  backplane  assumption.  The  maximum  errors  due  to  this  assumption  are 
of  the  order  of  40%. 


Table  B-2.  Comparison  of  the  effective  cross-coupling  resistances  from  a  floating 
backplane  setup  and  those  derived  from  a  grounded  backplane  for  a  multi-contact 
_ example  in  a  100pm  thick  substrate.  _ 


Rij 

From  floating 

Extracted  from  grounded 

Relative  Error 

backplane(fi) 

backplane(fi) 

Rl2 

225 

195 

13% 

Rl4 

333.33 

250 

25% 

R25 

409.83 

230 

43% 

R45 

322.8 

230 

28% 

As  the  bulk  thickness  reduces  below  200pm,  the  backplane  equipotential 
assumption  becomes  invalid.  For  processes  considered  in  this  work,  the  single  node 
approximation  of  the  backplane  is  valid  as  the  substrate  thickness  was  250pm. 
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APPENDIX  C 


The  choice  of  the  geometric  data  set  is  very  important  as  mentioned  in  Chapter  3.2.  The 
following  figure  represents  the  types  of  contacts  and  orientation  included  in  the  data  set. 


(a)  (b) 


(c) 


[=l  □  | 

(d)  (e)  (f) 


Figure  C-l.  Geometric  cases  and  orientations  considered  for  the  model  development. 

(a)  Identical  square  contacts,  (b)  Identical  rectangular  contacts  with  large  lengths,  (c) 
Identical  rectangular  contacts  with  large  widths,  (d)  Identical  rectangles  oriented 
differently,  (e),  (f)  Non  identical  contacts. 

The  geometric  test  cases  included  for  the  macromodel  development  are  listed  in  Table  C- 
1. 


Table  C-l.  Geometric  cases  for  the  model  development. 


Test  case  # 

W1 

LI 

W2 

L2 

1. 

1.2 

1.2 

1.2 

1.2 

2. 

2 

2 

2 

2 

3. 

3 

3 

3 

3 

4. 

4 

4 

4 

4 

5. 

5 

5 

5 

5 

6. 

6 

6 

6 

6 

7. 

8 

8 

8 

8 
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APPENDIX  D 

The  process  parameters  for  the  model  are  extracted  by  a  least-squares  curve  fitting 
procedure.  The  built-in  function  of  MATLAB  is  used.  This  section  will  explain  the 
procedure  for  extracting  the  process  parameters  for  the  GMD  model.  A  similar  procedure 
is  followed  for  obtaining  the  parameters  for  the  rest  of  the  model. 


The  objective  function  for  the  curve  fitting  procedure  is: 


min 


where  m  is  the  number  of  data  points,  Zmod  is  the  model  generated  Z  values  given  by  say 
Zmod  =  i  ( wi  +  w2 )  +  k2]e  k^gmd  anc)  zsim  js  the  EPIC  generated  Z  values.  The 


process  constants  for  the  model  are  extracted  from: 


k  =  arg  min 


where  k  is  the  vector  containing  the  process  constants  for  the  model  equation  and  arg  min 
returns  the  vector  k  that  minimizes  the  sum  of  squares  of  the  relative  errors.  The 
nonnalized  Z  values  are  used  for  the  curve-fitting  process  instead  of  the  simple 
difference,  to  make  the  procedure  more  robust. 
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APPENDIX  E 

The  model  development  procedure  used  115  contact  geometries.  However,  to 
extract  the  process  constants  of  the  macromodel  for  other  processes,  such  a  large  data  set 
may  not  be  available  and  hence  the  data  set  requirements  need  to  be  explored. 

From  the  available  data  set,  a  large  number  of  contact  geometries  with  small 
dimensions  (<20pm)  representing  different  orientations  are  excluded  to  create  a  new  data 
set  with  64  test  cases.  The  excluded  data  set  contains  contact  geometries  from  Test  case 
#14  to  #64  from  Table  C-l  (shaded).  The  process  parameters  for  the  model  equations  for 
the  new  data  set  are  extracted  using  the  least-squares  curve  fitting  procedure.  The 
comparison  between  the  old  process  constants  and  the  new  process  constants  are  shown 
in  Table  E-l. 


Table  E-l.  Comparison  of  process  constants  extracted  using  115  test  cases  and  64 

test  cases. 


Process  constant 

Value  with  115  Test  Cases 

Value  with  64  Test  Cases 

ki 

0.2322 

0.2384 

k2 

638 

637.34 

k3 

0.195 

0.195 

k4 

0.002 

0.002 

ks 

0.25 

0.25 

Oj 

9.5xl0"x 

9.5xl0"x 

(*2 

3xl0"6 

3.2xl0"6 

0.3 

0.5 

0.5 

a4 

4.6x1  O'4 

4.27xl0"4 
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It  can  be  observed  that  the  process  constants  do  not  undergo  significant  changes 
as  a  result  of  using  the  new  data  set.  Also,  the  model  equations  with  the  new  process 
constants,  were  validated  for  the  excluded  data  set  and  resulted  in  similar  error 
performance  with  the  original  model  equations.  Hence,  the  new  data  set  can  be 
considered  as  an  optimized  data  set.  The  contact  geometries  are  listed  in  Table  E-2. 


Table  E-2.  Geometric  test  cases  for  the  optimized  data  set. 


Test  Cases  # 

W1 

LI 

W2 

L2 

1. 

1.2 

1.2 

1.2 

1.2 

2. 

2 

2 

2 

2 

3. 

3 

3 

3 

3 

4. 

4 

4 

4 

4 

5. 

5 

5 

5 

5 

6. 

6 

6 

6 

6 

7. 

8 

8 

8 

8 

8. 

10 

10 

10 

10 

9. 

12 

12 

12 

12 

10. 

14 

14 

14 

14 

11. 

16 

16 

16 

16 

12. 

18 

18 

18 

18 

13. 

20 

20 

20 

20 

14. 

4 

2 

4 

2 

15. 

2 

2 

1 

2 

16. 

3 

5 

2 

3 

17. 

2 

5 

3 

4.5 

18. 

5 

10 

10 

4 

19. 

1 

30 

2 

2 

20. 

25 

25 

25 

25 

21. 

30 

30 

30 

30 

22. 

40 

40 

40 

40 

23. 

50 

50 

50 

50 

24. 

70 

70 

70 

70 

25. 

90 

90 

90 

90 

26. 

100 

100 

100 

100 

27. 

125 

125 

125 

125 

28. 

150 

150 

150 

150 

29. 

200 

200 

200 

200 

30. 

1 

30 

1 

30 

31. 

1 

40 

1 

40 

32. 

1 

100 

1 

100 

33. 

5 

40 

5 

40 

34. 

10 

40 

10 

40 
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APPENDIX  F 

In  order  to  compare  the  performance  of  the  macromodel  with  measurement  results 
the  entire  substrate  network  consisting  of  59  contacts  is  considered.  The  layout  of  the 
contact  geometries  is  shown  in  Figure  F-l.  Compared  to  the  layout  presented  in  Section 
5.4,  this  layout  includes  more  contacts. 
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Figure  F-l.  Layout  of  contacts  for  an  LNA  and  stepped  buffer  circuit  with  59  contacts. 


The  model  parameter  extraction  is  carried  out  using  the  optimized  data  set  and  the 
substrate  profile  shown  in  Figure  5.5  from  EPIC  simulations.  The  substrate  network  is 
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obtained  by  applying  the  paneling  approach  to  the  macromodel.  The  frequency  spectrum 
of  the  LNA  output  is  analyzed  to  estimate  the  substrate  noise.  The  model  is  reasonably 
accurate  when  compared  to  EPIC  and  measurement  results  as  shown  in  Figure  F-2. 


Figure  F-2.  Comparison  of  LNA  output  spectrum  using  macromodel,  EPIC  and 

measurements. 


